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Abstract 

The excitation spectrum of clean ferromagnet-superconductor bilayers is calculated 
within the framework of the self-consistent Bogoliubov-de Gennes theory. Because 
of the proximity effect, the superconductor induces a gap in the ferromagnet spec- 
trum, for thin ferromagnetic layers. The effect depends strongly on the exchange 
field in the ferromagnet. We find that as the thickness of the ferromagnetic layer in- 
creases, the gap disappears, and that its destruction arises from those quasiparticle 
excitations with wavevectors mainly along the interface. We discuss the influence 
that the interface quality and Fermi energy mismatch between the ferromagnet and 
superconductor have on the calculated energy gap. We also evaluate the density of 
states in the ferromagnet, and we find it in all cases consistent with the gap results. 

Key words: Energy gap, Proximity effect, Bilayers, Superconductivity 
PACS: 74.50.+r, 74.25.Fy, 74.80.Fp 



1 Introduction 



Through the mechanism of Andreev reflection [1], a normal (non-superconducting) 
material in contact with a superconductor acquires superconducting phase co- 
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herence between its particle and hole wave functions, even though the pair- 
ing interactions vanish in this normal material. The associated superconduct- 
ing correlations give rise to a finite value of the pair amplitude, F(r) = 
{^(r)ipi(r)}. This is known as the proximity effect. 

The proximity effect occurs also, in a modified way, if the non-superconductor 
is a ferromagnet. In this case, F(r) exhibits spatial oscillations. On the other 
hand, the pair potential itself, A(r) vanishes in the ferromagnet in the absence 
of attractive coupling. For conventional isotropic superconductors, A is spa- 
tially independent, and corresponds to the minimum excitation energy in the 
spectrum, or the energy gap, E g . For ferromagnet-superconductor systems, as 
in other inhomogeneous cases, A(r) depends on position, while the excitation 
gap for the whole system, E g , obviously does not. The proximity effect can 
therefore result in E g and A possessing a rather nontrivial relationship. In this 
paper we study this question for ferromagnet-superconductor (F/S) bilayers. 

Several works involving nonmagnetic normal metal-superconductor heterostruc- 
tures have discussed in various contexts, the presence or absence of an energy 
gap in the excitation spectrum. [2,3,4,5,6] For a normal metal layer of width 
greater than the coherence length in the superconductor, £ > but smaller than 
the phase coherence length L^, the appearance of a minigap of order of the 
Thouless energy has been reported, while if the normal layer is smaller than 
£oj the induced gap is of order A. It was found that while interface roughness 
plays a role, the minigap persistence depends chiefly upon interface quality. [7] 
These conclusions draw upon the results of calculations carried out within 
the quasiclassical formalism, which is applicable when certain conditions are 
satisfied. The main results are modified if the normal metal region is modeled 
as a non-Fermi liquid, where increased electron correlations may inhibit or 
completely destroy the minigap altogether. [8] 

For F/S structures, the magnetic exchange field in the ferromagnet results in 
quasiparticle spin splitting, and consequently the pair amplitude rapidly de- 
cays and oscillates over a characteristic length scale £p, in which electron-hole 
pairs remain coherent. This modulation in phase of the pair amplitude has 
been shown to induce corresponding oscillations in the local density of states 
(DOS). [9,10,11,12] Impurities in the magnet were shown to play a significant 
role in the determination of thermodynamic quantities. [9,11] For ballistic F/S 
structures, the tunneling DOS for clean ferromagnet layers has been calculated 
for a wide range of exchange fields, [12] and the local DOS has been found to 
exhibit a gapless state for all positions in the magnet. [10,13] Configurations 
consisting of weak or thin ferromagnetic layers where dp/t^F <C 1, m &y lead 
to significantly different results, [14] as we shall see below. The finite geometry 
can result in the quasiparticle amplitudes undergoing coherent Andreev and 
normal reflection (at the ferromagnet-insulator interface) in such a way that 
the minimum value of the excitation spectrum is nonzero due to the relatively 
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slow characteristic decay of the pair amplitude. Thus, a gap develops. Associ- 
ated with this, when the dimensions of the ferromagnet are reduced to a scale 
comparable to £p, the greater confinement of the pairing correlations in the 
ferromagnet restricts the pair amplitude so that it no longer changes sign, and 
the DOS does not exhibit the corresponding oscillations. 

It is generally accepted that the quasiclassical and dirty limit equations are 
satisfactory for studying many aspects of inhomogenous superconductivity. 
If however, the ferromagnet has dimensions of order of or smaller than the 
elastic mean free path £, so that quasiparticles undergo very few scattering 
events, the problem is more suitably solved within the ballistic regime. Fur- 
thermore, the inherent geometrical effects associated with ferromagnet layers 
only a few atomic spacings thick require a full microscopic theory that does not 
coarse-grain over length scales of order of the Fermi wavelength X F . Another 
constraint imposed upon the quasiclassical approximation is exhibited by the 
Andreev[l,15] equations, whereby values of the quasiparticle momentum par- 
allel to the interface, and comparable to the Fermi momentum, are crudely 
approximated. We will see below that in fact these quasiparticle states play a 
prominent role in the determination of the gap threshold for F/S systems. 

Thus, we examine here the electronic spectrum of F/S bilayers and investigate 
under what conditions an energy gap can exist in such structures. The prob- 
lem will be solved within the framework of the self- consistent solutions of the 
Bogoliubov de-Gennes[16] (BdG) equations. We will calculate self-consistently 
the DOS, and excitation spectra of a thin ferromagnet film adjoining a bulk su- 
perconductor. This microscopic approach complements existing quasiclassical 
theories and serves to round the current knowledge. 

In the next section, we introduce the method and geometry used in our model. 
In Sec. 3, the influence of the interface transparency, the magnitude of the 
exchange field, and the mismatch in Fermi levels on the energy gap and local 
DOS, will be presented. 



2 Geometry and model 

We consider a heterojunction of total length d in the z-direction, consisting 
of a ferromagnet of thickness dp and a superconductor of thickness ds- The 
planar interface is at z = d F , and the free surfaces at z = and z = d are 
specularly reflecting. Upon taking into account the translational invariance 
in the semi-infinite x — y plane, one can immediately write down the BdG 
equations[16,10] for the spin-up and spin-down quasiparticle and quasihole 
wave functions (u^v^), 
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+ e ± + W(z) - E F {z) vi(z) + A(z)ul(z) = e n vi(z), 



+ e ± + W(z) - E F (z) ul(z) + A(z)v l n (z) 




(1) 



(2) 



where e± is the transverse kinetic energy, e n are the quasiparticle energy eigen- 
values, h (z) = h Q(z — dp) is the magnetic exchange energy, and A(z) is 
the pair potential. Scattering at the interface is modeled by the potential 
W(z) = W5(z — (If), where W is the barrier strength parameter. We take the 
quantity Ep(z) to equal Efm in the magnetic side, so that Epi = Epu + ho, 
and Epi = Epm — ho, while in the superconducting side, Ep(z) = Ep. From 
the symmetry of the problem, the solutions for the other set of wavefunctions 
(u^vl) are easily obtained from those of Eqns. (1,2) by allowing for both 
positive and negative energies. The BdG equations are supplemented by the 
self consistency condition for the pair potential, 



where T is the temperature, g(z) is the effective coupling constant, describing 
the electron-electron interaction, and ojd is the Debye energy. 



3 Results 

In the following, we consider a thin ferro magnet layer adjacent to a bulk limit 
(ds ^> £o) superconductor. To focus on the parameter range which may yield 
a gap in the energy spectrum, the thickness of the ferromagnet considered 
will not exceed the coherence length in the superconductor, i.e., dp < £ - We 
introduce the dimensionless parameter Z to characterize the barrier strength, 
with Z = mW/kp, and the measure of mismatch between the Fermi levels is 
given by the ratio A = EpujEp. The requirement of self-consistency in the 
problem excludes the possibility of an analytical solution. The BdG equations 
(1,2) are therefore solved numerically using the approach [10, 17] outlined in 
the Appendix. We assume the low temperature regime (we take T/T c = 0.02), 
and a superconducting coherence length £ of kp£, = 50, where kp is the Fermi 
wavevector of the superconductor. 

The energy gap, E g , is the minimum binding energy of a Cooper pair, and 
its existence will influence various thermodynamic measurements, including 
heat capacity and thermal conductivity. For other configurations, such as long 
superconductor-normal metal-superconductor junctions, the gap in the metal 
and the bulk gap in the superconductor can further inhibit thermal transport. 
The process of finding E g is computationally demanding, as it involves cal- 




]T [ul(z)vi(z) + ui(z)vl(z)} tanh(6 n /2T), 



(3) 
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Fig. 1. Variation of the normalized energy gap in the ferromagnet with the dimen- 
sionless ferromagnet thickness kpd.F- The coherence length is given by &f£o = 50. 
The different curves correspond to the four values of the magnetic exchange energy 
labeled in the figure. 



culating self-consistently the entire eigenvalue spectrum and then finding its 
minimum, for each value of d F - In Fig.l we illustrate our results for the de- 
pendence of Eg on the thickness of the ferromagnet layer. These results are for 
Z = and A = 1. Results for four different exchange fields in the ferromagnet 
are presented. The top dashed curve corresponds to ho/A Q = 0, while the sub- 
sequent lower curves are for finite values of the exchange field as labeled in the 
figure. All curves start at E g = A for d F = 0, corresponding to the expected 
result for a single superconductor in the bulk limit. As the ferromagnet layer 
grows, Eg decays monotonically towards a gapless superconducting state, at a 
rate that depends strongly on the strength of ho- In the limit of zero exchange 
field (normal metal), the gap ceases to exist when dp ~ £ - This thickness at 
which Eg is destroyed is much smaller than the proximity length character- 
izing the usual very slow decay of the pair amplitude in the metal. [18] This 
trend continues for finite ho, with E g vanishing much more rapidly than the 
characteristic decay length of the pair amplitude in the magnet. To illustrate 
this, recall[17,19] the approximate expression for the characteristic length of 
decay of the pair amplitude in the ferromagnet: £p ~ {k F ^ — kpi) -1 , which 
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Fig. 2. Part of the normalized self consistent quasiparticle spectrum as obtained 
for given values of the normalized transverse energy e±/Ep. The parameter /io/Ao 
is chosen to be 0.8, and IzfcIf = 15. The region | e n / Ao | > 1 is dominated by the 
continuum sates, these are suppressed for clarity. 

can be rewritten as, 



Thus for example, the h /A = 4.0 curve in Fig. 1 shows the energy gap 
vanishing at /c^m^f ~ 7.5, while the pair amplitude decays over the larger 
length scale /cfm£f ~ 20. The discrepancy is even larger if one recalls that 
the expected decay length might be thought to be a factor of 2tc larger. For 
the largest field shown (h /A = 8), the superconducting correlations are 
further inhibited, failing to generate a gap in the energy spectrum except for 
ferromagnet nanostructures satisfying /cpd^ < 5. 

To provide additional perspective on the spectral features of the proximity 
effect, we show in Fig. 2 a relevant portion (|e n /A | < 1) of the calculated 
eigenvalues as obtained (under the same conditions as in the previous Figure) 
for given transverse energies, e±. For a homogeneous superconductor, a gap 
separates the continuum states that exist for |e n /Ao| > 1. With the inclusion 
of a ferromagnet, in addition to the scattering states there exists a discrete 
quasiparticle excitation spectrum in the range |e n /A | < 1 (see Fig. 2). The 
bound states are visibly asymmetric about the Fermi level due to the effects 
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the exchange field has on the particle-hole amplitudes. The existence of a 
finite number of bound-state branches of localized quasiparticles is also char- 
acteristic of antisymmetric domain walls. [20] In Fig. 2 we can further see that 
the minimum in |e n |, which corresponds to the energy gap seen in Fig. 1 for 
kFd-F = 15, where E g pa 0.09A , is sharp, and that it occurs for quasiparticles 
whose in-plane momentum is close to the Fermi level (e±/Ep — 1). This is due 
in part to those quasiparticles not coupling to those states responsible for su- 
perconductivity. This result is typical of what we find for other exchange fields, 
the only difference being the value of E g [see Fig. 1]. It is therefore those excita- 
tions with a dominant momentum component parallel to the interface, that are 
significant in contributing to the filling in of the gap. In the Andreev or quasi- 
classical approximation scheme, these quasiparticles have "trajectories" which 
are not treated accurately, and those with e±/E F > 1 are neglected altogether. 
This discrepancy has also been revealed [15] in normal metal-superconductor 
systems, where normal reflection dominates the Andreev reflection process for 
quasiparticles with large momentum parallel to the interface. 

The previous results in Fig. 1 showed the strong influence the exchange energy 
has on E g . Another quantity of great experimental interest is the local density 
of one particle excitations in the magnet. Current experimental tools such as 
the scanning tunneling microscope (STM) have atomic scale resolution, and 
make this quantity experimentally accessible. When well defined quasiparticles 
exist, the tunneling current is simply expressed as a convolution of the one- 
particle spectral function of the STM tip with the spectral function for the F/S 
system. [21] The resultant tunneling conductance, which is proportional to the 
density of states (DOS), is then given as a sum of the individual contributions 
to the DOS from each spin channel. 

N(z,e) = N 1 (z,e) + N l (z,e), (5) 
where the local DOS for each spin state is given by 

Nfa *) = - H{Wn{z)] 2 f{e - e n ) + \vl{z)f f'(e + e n )}, (6) 

n 

Ni(z, e) = - £{kW] Y( e - e n ) + [vi(z)ff\e + e n )}. (7) 

n 

Here thermal broadening is accounted for by the derivative of the Fermi func- 
tion, /'(e) = df/de. To investigate further how the previous results interrelate 
with the local DOS (Eq. 5), we show in Fig. 3 the local DOS N(e), defined 
as N(z, e) integrated over the ferromagnetic layer. This quantity is shown for 
two exchange fields, differing by a factor of ten. The ferromagnet has a thick- 
ness of d,F = 9kp 1 : from Fig. 1 this corresponds to an energy gap of E g pa 0.5 
for h /A = 0.8, and an absence of an energy gap for h /A = 8. This is 
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Fig. 3. Density of states in the ferromagnet (see text) for two exchange energies 
and a fixed width /^dp = 9. The bold curve corresponds to Jiq/Aq = 0.8, and 
the dashed curve to h /A = 8. The DOS is normalized to N(0), the DOS of the 
superconductor in the normal state, and the energies are normalized to Ao- The 
inset reveals the same quantities, magnified to clarify the gap region. As in previous 
Figures, we assume no interfacial scattering and no Fermi wavevector mismatch: 
Z = 0, A = 1. 

consistent with the DOS calculations, as Fig. 3 illustrates: for the smaller ex- 
change field, there is a complete absence of states for e/A < 0.5, while the 
DOS is gapless at all energies for h /A = 8. The observed energy gap for 
ho/Ao = 0.8 is a result of the underlying competition between quantum size 
effects and constructive interference between the spin-split quasiparticle am- 
plitudes. As expected, the particle-hole asymmetry is more prevalent in the 
stronger magnet (within the subgap region), as the two uppermost peaks are 
clearly shifted relative to one another. 

Next we consider the effect that an interface barrier or Fermi energy mismatch 
has on the electronic structure. Experimentally, some degree of Fermi energy 
mismatch is unavoidable, particularly with alloys, and the transparency of 
the interface can be degraded due to a thin oxide layer. In Fig. 4 we illustrate 
the normalized energy gap E g as a function of kpMdF (where kpu is the 
wavevector corresponding to Efm), f° r two A and Z values that each differ 
by a factor of three. We fix the exchange field at ho = 0.8A , keeping in mind 
that higher fields obey similar trends so that generality is not lost. Also, in 
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Fig. 4. Influence of interface scattering and Fermi energy mismatch on the energy 
gap as a function of the ferromagnet thickness. Both the effects of increasing the 
mismatch (decreasing A from unity, see text) and increasing scattering strength 
(increasing Z) serve to destroy the energy gap E g at an increased rate. 



order to isolate each of the effects, whenever A is different from unity, Z is 
zero, and vice versa. The figure shows that increased mismatch (smaller A) or 
increased interface scattering (larger Z), induces the same trend of reducing 
Eg as a function of hpM^F ■ The functional form of the graphs are distinct 
however: for A < 1, the tail-end decay is slower than the near linear decay in 
E g exhibited in the case of finite Z. Further differences in the effects of the 
two parameters are even more evident in the DOS spectra. To this end, Fig. 5 
shows the DOS in the ferromagnet for two different values of the mismatch 
parameter A, and one value of Z. Upon comparison with the previous Fig. 4, 
we again find consistency between the DOS and E g calculations: the energy 
gap is reduced at kpMdF = 10 for both A = 0.2, and Z = 0.2. This is reflected 
in the DOS, where Fig. 5 illustrates a small but finite number of states centered 
about the Fermi level. Although the quasiparticle spectrum has just turned 
gapless for A = 0.2, and Z = 0.2, the DOS signature is substantially different, 
the most prominent feature being the larger number of subgap bound states 
for A = 0.2, which can be attributed to the larger Fermi energy mismatch, 
and subsequently, increased normal reflections. 
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Fig. 5. Local DOS in the ferromagnet of width kpudF = 10. The bold curve corre- 
sponds to Z = 0.2, and the dotted and dashed curves are for the cases A = 0.6, and 
A = 0.2 respectively. For the given dimensions of the magnet, the inset reveals that 
only when A = 0.6 does a gap survive in the energy spectrum. At higher energies 
the curves are clearly displaced from one another by a factor of y/X, due to the 
normalization. 



4 Conclusions 



We have self-consistently calculated the electronic structure for a F/S het- 
erostructure consisting of a thin ferromagnet layer adjoining a bulk supercon- 
ductor. Our fully microscopic method revealed the dependence of the transi- 
tion to a gapless superconducting state on the relevant physical parameters 
such as exchange field energy, interface transparency, and Fermi energy mis- 
match. The energy gap decreases with the width of the F layer at a rate that 
depends upon the dimensionless ferromagnet width hpu^F- The characteristic 
length for gap suppression is much smaller than the length scale describing the 
decay of the pair amplitude in the ferromagnet: thus E g is suppressed more 
rapidly than the pair amplitude. Fermi level mismatch and interface barrier 
both serve to reduce the gap further, however the functional forms of the re- 
spective reductions are different. The local DOS in the ferromagnet is in all 
cases consistent with the above observations. The self-consistent quasiparticle 
spectra reveals that the gap onset is due to quasiparticle states with a large 
momentum component parallel to the interface. These states however, occupy 
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a relatively small region of momentum space, and when calculating integral 
quantities over the full solid angle, good agreement between our exact theory 
and quasiclassical theory might be possible. 

We have focused here on the spectral properties of the ferromagnet, in the 
clean limit. This regime is most appropriate for layers whose dimensions do 
not exceed the mean free path (d F < £). Since, we have considered only thin 
layers, this assumption is consistent with our calculations. For situations where 
the effects of impurities or inelastic effects may be important, a finite DOS 
at low energies would likely arise. The general method used here however has 
found good overall agreement with experiments involving thin nonmagnetic 
normal metal-superconductor bilayers, and bulk F/S structures. [17] Further- 
more, previous calculations[13,22] within the ballistic limit have also found 
agreement with experiment. Other complexities may arise in actual experimen- 
tal conditions, such as grain boundaries and domain walls not accounted for 
here. The fabrication of relatively clean heterostructures, [23,24] and structures 
with highly transparent interfaces [25] is possible however. Our calculations are 
therefore realistic, and the results should be reproducible experimentally. 



A Numerical method 



We solve Eqns. (1,2) by expanding the quasiparticle amplitudes in terms of 
a finite subset of a set of orthonormal basis vectors, u^z) = Y,q u nq < t>q(. z )i 
and v^{z) = *52q v nq < t>q( z )- We use the complete set of eigenf unctions 4> q {z) = 

(z\q) = ^2/dsm(k q z), where k q = q/nd, and q is a positive integer. The finite 
range of the pairing interaction ujd permits the number N of such basis vectors 
to be cut off in the usual way. [10] Once this is done, we arrive at the following 
2N x 2N matrix eigensystem, 



H+ D 
D H 



(A.l) 



where ^ = (u nl , . . . , u nN , v^\, • • • , v^ N ). The matrix elements H^ q , connecting 
<p q to (p q i are constructed from the term found in brackets in Eq.(l), 



U qq' 



k^ 

2m 



1 d 2 



2m dz 2 



+ e ± + W(z) - E F {z) 



+ e± 



d dp 

5qq> + J dz(f) q (z)W(z)(/) q/ (z) - E F] J dz (f) q (z)<p ql (z) 
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u 

-E F J dz <p q (z)(j) q/ (z). 

dp 



(A.2) 



The expression for H , is calculated similarly using Eq.(2). The "off-diagonal" 
matrix elements D qq i are given as, 



a 

D qq , = {q\A(z)\q'} = J dz c/ )q (z)A(z)(f> q '(z). 

dp 



(A.3) 



After performing the integrations, Eq.(A.2) can be expressed as 



2W 



H q , = —— sm(k q d F ) sm(k q id F ) 



E 



Ft 



d 

E F 

~d~ 



d 



sia[(kgi — k q )d F ] sin[(A; g / + k q )d F ] 



\kqi k q ) 



sin[(A; 9 / + k q )d F ] sm[(k q > — k q )d F ] 

(kg' + kg) (k q ' - kg) 

while the diagonal matrix elements are written as, 



(A.4) 



r kl 



n qq> 



d s + 



2m 
E F 



W 



1 — cos(2k q dp) 
sm(2k q dp) 



E 



FT 



d 



di 



sin(2k q dp 
2k a 



2k a 



q = q 



The self-consistency condition, Eq.(3) is now transformed into, 
ttA 



r p,p' i 



(A.5) 



11^ 7?^ I 7/^ 77^ 
Uj np u np' ' Uj np u np' 



sin(kpz) sin(kp'z) tanh(e„/2T), (A. 6) 



where A = g(z)N(0), and iV(0) is the DOS for both spins of the superconduc- 
tor in the normal state. The quantum numbers n encompass the continuous 
transverse energy e±, and the quantized longitudinal momentum index q. The 
matrix eigensystem Eq. (A.l) and the self-consistency condition (A. 6) are 
then solved numerically, using an iterative scheme developed and described in 
earlier work[10,17]. 
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